Compression as a tool to detect Bose glass in cold atoms experiments 
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We suggest that measuring the variation of the radius of an atomic cloud when the harmonic 
trap confinement is varied makes it possible to monitor the disappearance of the insulating Mott 
phase of an ultracold atomic gas trapped in a disordered optical lattice. This paves the way for an 
unambiguous identification of a Bose glass phase in the system. 

PACS numbers: 03.75.Lm,67.85.Hj,03.75.Kk 



Condensed matter models have found recently a won- 
derful testbed in cold atoms in optical lattices physics 0, 
■ Cold atoms allow for an optimal control of parame- 
ters of the system: by changing intensities or detunings 
of laser beams, one may modify the depth of the opti- 
cal lattices; interactions between atoms can be controlled 
via external magnetic field and Fcshbach resonances. An 
example is the Bose-Hubbard tight binding model 
which can be realized using ultracold atoms and where 
the quantum phase transition between a supcrfluid (SF) 
phase and a Mott insulator (MI) phase can be observed 
by varying the depth of the optical lattice as predicted 
in a seminal work 0] and realized few years later 

Controllability of the system makes studies of disor- 
der induced effects particularly interesting. Disorder (or 
pseudo-disorder) may be created in cold atoms systems 
in a rcpeatable way using optical potentials created by 
laser speckles or multichromatic lattices [H, 0] ■ This ex- 
citing possibility attracted soon a lot of research [f| - 
Anderson localization being one of the main targets for 
weakly or non- interacting atoms [H, [13] • Even more ex- 
citing is another regime of strongly interacting bosons 
in a disordered potential. Seminal studies of the disor- 
dered Bose-Hubbard model 0, EH revealed the existence 
of a novel insulating phase called Bose glass (BG) phase. 
Contrary to the MI, the BG is characterized by gapless 
excitation spectrum and is compressible. As far as we 
know it has not been observed in a "traditional" con- 
densed matter settings. Its unambiguous observation is 
certainly an important milestone yet to be seen. 

In the first attempt to produce a BG with ultracold 
atoms, a bichromatic quasi-disordered optical lattice was 
used |12|. The authors measured both the absorption 
spectrum of the system and the long-range spatial coher- 
ence. They observed a smearing of the absorption peaks 
in the presence of "disorder" as well as a decreased long- 
range coherence, which they interpreted as a manifesta- 
tion of the existence of a BG phase. The interpretation 
of absorption spectra is however difficult: as discussed 
by us elsewhere [HI , the initial state was not the ground 
state of the system (this would be the case if the lattice 



were ramped up adiabatically as supposed in Ref. [13]). 
Also the amplitude of the lattice modulation was strong, 
so that the absorption was not in the linear regime. 

Another very recent attempt [IH has shown - in the 
three-dimensional case - that the presence of the disorder 
leads to a significant decrease of the condensate fraction 
both for SF or coexisting SF and MI phases. This mea- 
surement follows in fact the original proposition of Q to 
address the disappearance of the condensate fraction as 
a possible signature of the BG presence. Here, again the 
preparation of the initial state is a key question. 

While both these experiments are important studies of 
strongly interacting bosonic systems in the presence of 
disorder, it is desirable to have a clear signature of a BG 
phase. Of special interest is a direct measurement of the 
compressible or incompressible nature of the system. An 
interesting possibility is to access the central density of 
the atomic cloud and measure the changes of that den- 
sity when the external t rapp ing potential is varied (trap 
squeezing spectroscopy) [15| . That measurement supple- 
mented with coherent fraction measurement provides an 
interesting and clear proposition for an "ideal" experi- 
mental procedure. There are two requirements. One is 
to make the trap geometry independent from laser beams 
forming the lattice, so that one can vary the trap fre- 
quency independently of the lattice height. The second 
one is to use an additional focused laser to monitor the 
density in the center of the trap. The latter while feasible 
seems quite difficult. The aim of this letter is to propose 
a much simpler experimental scenario, which has already 
been used in demonstrating the Mott phase for a cloud of 
fermions [lij]. We propose to measure the radius of the 
atomic cloud while changing the trap frequency. This 
provides a simple and direct measurement of the com- 
pressibility of the system! More precisely, it makes it 
possible to monitor the appearance and disappearance 
of incompressible phases. Complemented with measure- 
ments of the long-range phase coherence, it would allow 
for an unambiguous characterization of the MI-BG-SF 
phase diagram. 

This method works beautifully for fermions fl6j | where 
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in the Mott state one has at most one fermion per spin 
state. In such a case the radius of the cloud becomes 
practically independent of the trap frequency, clearly 
demonstrating the incompressibility of the MI state. The 
situation is quite different for bosons where in typical ex- 
periments 0, [13 one may have up to three bosons per 
site of the optical lattice and the density profile in the 
trap resembles that of a wedding cake - see e.g. fl7| . 
Between Mott regions with single, double, or triple oc- 
cupancies, there are supcrfluid "shells" . Those lead to a 
finite compressibility of the total sample and could make 
radius measurements useless. As shown below, this is not 
a major problem. 

Let us start with a one-dimensional Bose-Hubbard 
model in the presence of a trapping potential and a di- 
agonal disorder. The Hamiltonian is 



H = - J £ b]b f 



U 



fij {hj - 1) + ^ , 



(1) 



where bj is the destruction (creation) operator of one 

particle in the j-th site, fij = b^bj is the number operator, 
and (j, f) indicates the sum on nearest neighbors. U is 
the interaction energy and J the hopping energy. The 
energies at sites, ej are given by the sum of the energy 
shift due to the the harmonic potential and the additional 
disorder: 



ej = ^7iu) 2 a 2 (j - j ) 2 + XjUA 



(2) 



where m is the particle mass, a the lattice spacing, cu the 
trap frequency and jo the position of the trap center. A 
is a dimensionless parameter measuring the strength of 
the disorder (in units of the interaction energy), while Xj 
is a (pseudo-) random variable. We consider two different 
types of disorder. For a truly random disorder, the Xj's 
are chosen as independent variables with uniform dis- 
tribution in the interval [-1,1]- For a secondary optical 
i is a sine function with incom- 



lattice as used in [12J, 
mensurate frequency: 
random correlated variables with distribution 1 



P(x) = 



sin(Aj) resulting in pseudo 



(3) 
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The system properties in such a pseudo-random disorder 
may differ from the truly random situation 2, 2^, 21 1 . 

We employ the parameters originating from the Flo- 
rence [l2j experiment with the exception that we assume 
the possibility of independent change of trap frequencies. 
In particular we concentrate on the deep Mott regime 
(J/U « 0.027). 

As shown below, in this regime - and as long as we 
look at the compressibility of the system and not at 
long range phase coherence - we can use a local de nsity 
approximation (LDA) to describe the system [l6|, (Hf- 



1.2 


1 

0.8- 
0.6 
0.4, 




2 4 6 



0.5 



1 1.5 



2 2.5 3 

LO 



Figure 1: Scaled r.m.s. radius of the one-dimensional atomic 
cloud as a function of the scaled trap frequency for a random 
uniform disorder with A=0.2. The kink around loq — 1 is a 
signature of the appearance of a new phase. The solid line is 
obtained using the local density approximation (LDA). The 
circles are the results of a quantum calculation with the TEBD 
algorithm for 61 particles and J/[/=0.027. The inset shows 
the scaled r.m.s. radius on a wider scale, in the absence of 
disorder. One clearly sees the additional kinks related to the 
appearance of doubly and triply occupied MI phases. 



This amounts at neglecting tunneling between neighbor- 
ing sites and assuming a Fock state at each site with the 
occupation determined by the local chemical potential 
fi — ej. Determining the shape of the atomic cloud is a 
simple minimization procedure for the total energy, i.e., 
the sum of independent contributions for various sites, 
constrained by a fixed total number of particles. 

In the absence of disorder and at low trapping fre- 
quency, the ground state is a pure MI with an unit occu- 
pation number at all sites around the trap center. If we 
define the r.m.s radius of the trap (in units of the lattice 
spacing a) as: 



R=^(ri)-(r)2 = 



N 



N 



(4) 



where N is the total number of particles, it is clear that, 
for large N and low trap frequency, the radius will be 
R c = N/2y3 independently of the trap frequency, a 
clear-cut manifestation of the incompressibility of the MI 
phase. When the trap frequency is increased, the energy 
of the outest particles m(jj 2 a 2 N 2 /8 increases until the 
point where it is cheaper to pay the interaction energy 
U and transfer the particle at the trap center, creating 
doubly occupied sites. Beyond this point, R decreases 
with frequency, implying global compressibility. In the 
large N limit, a straightforward calculation shows that 
the parameters obey scaling laws. The critical frequency 

is (from the two estimates above) uj c = J -^rjj ■ We thus 
define the scaled frequency: ojq = uj/lu c and the scaled 
r.m.s. radius, Rq = R/R c . 

A plot of Rq versus ujq is shown in the inset of Fig.[l] for 
iV=151 in the absence of disorder. For shallow trap, the 
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radius shows a pronounced plateau. In this range of fre- 
quencies, all particles are in the Mott phase with an unit 
filling - the plateau is a direct manifestation of incom- 
pressibility. This regime resembles most the fermionic 
case fl6| . At u>q=1, a sharp kink indicates the moment 
when the sample becomes compressible - at this point a 
double occupancy appears at the centre of the trap (this 
can be visualized looking directly at the occupation of 
sites). The compressibility becomes smaller for larger 
frequency reaching a second kink at dq = 1 + v2 ~ 2.414 
- here a triple occupancy in the center of the trap is 
born. A careful inspection reveals even the third kink at 
u = l + \/2 + \/3~ 4.146. From this plot, it is clear 
that measuring the global quantity R vs. trap frequency 
makes it possible to monitor the appearance of the suc- 
cessive MI pha ses. Let us note that for the Florence 
experiment [12J the effective harmonic binding (coming 
from the trap and lasers' transverse profiles) is 75 Hz, 
corresponding to u n sa 3.44. 

In the presence of a random disorder, the plateau dis- 
appears, showing a non zero compressibility of the sam- 
ple, as seen in the main panel of Fig.[TJ At low frequency, 
this is due to the appearance of a compressible BG phase 
in the external part of the cloud (where the average oc- 
cupation number is between and 1). The kink around 
luq=1 is, however, a robust feature. It is due to the for- 
mation of first a BG with occupation number between 1 
and 2 at the trap center, followed by the birth of a n=2 
MI phase. The kink itself is the signature of the appear- 
ance of new phases. The residual tunneling between the 
neighbouring sites has almost no effect. To prove this 
statement, we use the TEBD algorithm [22| (also known 
as t-DMRG algorithm j23|) with imaginary time prop- 
agation to produce the quasi-exact ground state of the 
Bose-Hubbard Hamiltonian in the presence of disorder 
from which the corresponding r.m.s. radius easily com- 
puted. The result, averaged over several realizations of 
the disorder, is shown in Fig. [TJ It is almost indistin- 
guishible from the result of the LDA approximation. We 
have checked that, for other disorder strengths - but still 
in the low tunneling regime - there is a similar agreement 
between the exact TEBD result and the LDA approxi- 
mation, which we will use in the following of this paper. 
Most probably, this is because tunneling creates coher- 
ence between neighbouring sites and smoothes the steps 
of the wedding cake, but does not induce macroscopic 
transfer of particles over long distances and thus only 
marginally affects the r.m.s. radius. From the point of 
view of compressibility, whether a BG or a SF phase is 
formed makes little difference. 

The MI phase is expected to disappear at A=0.5 in 
the limit of small tunneling @. In Fig. 1 we show 
Ro vs. ojq for increasing disorder strengths. The 
kink is very clearly visible even at values very close to 
A=0.5, showing that the proposed simple method al- 
lows one to monitor appearance and disappearance of 
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Figure 2: (Color online) Scaled r.m.s. radius of the one- 
dimensional atomic cloud as a function of the scaled trap fre- 
quency for a random uniform disorder. The successive curves 
(from left to right) are for an increasing disorder strength 
A = 0, 0.05, 0.2, 0.4 and 0.5. The kink in the curves disap- 
pears for the limiting case of A=0.5 (dashed line) beyond 
which the MI phase no longer exists. The inset shows the 
dimensionless global compressibility, eq. (|5j, which displays a 
marked minimum related to the existence of the incompress- 
ible MI phase. 

various phases. At the critical disorder A=0.5, the 
MI phases completely disappears and a smooth curve 
is obtained. It can be shown [3] that its equation is: 

R = 2- 1 /3 3 5/6 5 -l/2 CJ -2/3^ 

The kink is even more visible if one considers the com- 
pressibility of the system, that is the derivative of the 
radius with trap frequency. In [l6j | . a global compress- 
ibility k r ; = - ^J 2^2/ 2 ) is defined, where the powers 
of R are chosen to obtain well-defined quantities in the 
thermodynamic limit N — » oo. As it has the dimension of 
the inverse of an energy we prefer to use a dimensionless 
global compressibility K by multiplying by U. It has 
a simple expression in terms of scaled quantities: 

K= U OR 3_dRo 

R 3 d{muo 2 a 2 /2) 2iJgw du 1 ' 

It is shown in the inset of Fig. [5] for various disorder 
strengths. While it vanishes below u->o=l for zero disor- 
der, it displays a well marked minimum for < A < 0.5, 
directly related to the existence of an incompressible MI 
phase. At A = 0.5, the MI phase and the minimum dis- 
appear. Note that the observed behaviour is quite similar 
to the one observed for fermions in [lj|, the role of tun- 
neling being here replaced by disorder. 

A similar behaviour is observed for the three dimen- 
sional isotropic Bose-Hubbard model. The results arc 
shown in Fig. [3] The singularity is slightly different 
(the compressibility has no discontinuity, but a very fast 
increase), but the essential property: the existence of 
a kink in the Ro vs. luq curve, or equivalently, the 
marked minimum in the compressibility, is present as 
well. The scaling laws are of course different in three di- 
mensions: we now have R c = 3 5 / 6 5~ 1 / 2 2 _1 / 3 7r~ 1 / 3 A 1 / 3 
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Figure 3: (Color online) Same as Fig. [2] for the three- 
dimensional disordered isotropic Bose-Hubbard model. The 
behaviour is similar, with a kink in the curve around uoa=l 
and a marked minimum of the compressibility (inset). 




Figure 4: (Color online) Same as Fig. [2] but for the pseudo- 
random "disorder" created by a secondary incommensurate 
lattice: A=0 (solid line), A=0.2 (dotted line), A=0.4 (dashed 
line), A=0.5 (dash-dotted line beyond which the kink and the 
MI phase disappear). 



and u c = 2 7 /6 7r i/3 3 -i/3 iV -i/3^ [/ / ma 2 and the 3/2 co- 
efficient in eq. ([5]) must be replaced by 5/6. The appear- 
ance of the second MI phase is clearly visible as a dip in 
the compressibility around wo = 1.5. 

Consider now the pseudo-disorder, eq. employed in 
the Florence experiment [l2T ] . We have performed similar 
calculations, averaging over the relative phase of the two 
lattices. An overall similar behaviour, shown in Fig. [4) is 
observed, a pronounced kink yields the frequency when 
the double occupancy emerges at the centre of the trap. 
With increasing disorder, the plateau shrinks, becomes 
tilted, the kink moves towards smaller frequencies and 
eventually disappears at A=0.5. This translates in the 
second lattice depth equal to S2 ~ 1.014 (in the notations 
of [13]), which nicely matches the estimate for the disap- 
pearance of the Mott phase and the gap in the absorption 
spectrum (l2| . 

Note that another singularity (sharp peak in the com- 



pressibility) is visible at lower frequency and large dis- 
order. This is due to the peculiarity of the distribution, 
eq. ((3j) , diverging at x = ±1. Those, however, may be 
easily identified and distinguished from the "main" Mott- 
end kink. 

It has been argued that the MI regions may be strongly 
affected by temperature effects [25| . Our calculations are 
limited to T — 0. Study of the effects of the temperature 
as well as of large tunneling on the phenomenon discussed 
here are in progress. 

In summary, we have shown that a simple measure- 
ment of the radius of the atomic cloud may provide a clear 
identification of the disappearance of the Mott phase 
and, in the presence of disorder, may help to unravel the 
BG presence. The disappearance of the MI is correlated 
with vanishing of the kink in the radius-frequency de- 
pendence. DD acknowledges support by IFRAF, JZ ac- 
knowledges support by Polish Ministry of Education and 
Sports (2008-2011). This work is realized within Marie 
Curie TOK scheme COCOS (MTKD-CT-2004-517186). 
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